La correcta asimilación de observaciones de satélite, en particular radianzas, depende en gran medida de la estimación y corrección del bias asociado a cada observación. Si el bias no es corregido podría afectar el análisis y los pronósticos generados generando impactos negativos.

Uno de los desafíos en la corrección de bias para asimilación regional es la cantidad de observaciones disponibles en cada ciclo de asimilación lo que depende a su vez de la ubicación y orbitas de los distintos satélites. En la región de Sudamérica los distintos satélites polares pasan una o dos veces por día y el dominio escaneado es variable. A esta restricción se suma la disminución de observaciones disponibles debido al límite superior de los modelos regionales.

El sistema GSI-LETKF incorpora un metodo de corrección de bias que estima y actualiza coeficientes \(\beta_i\) que se aplican a una regresón lineal previo a la estimación del estado final o análisis (Miyoshi, Sato, and Kadowaki 2010). Esta estimación se calcula aplicando la ecuación:

\[ \delta\beta = (B_{\beta}^{-1} + PR^{-1})^{-1}PR^{-1}[y-H(x)-P^{T}\beta]\]

Esta formulación es equivalente a la utilizada por la versión variacional de GSI y análoga a la ecuación del cálculo del análisis con la inclusión de la corrección de bias. En la ecuación \(\delta\beta\) se asume diagonal y sus elementos se inicializan con el valor 0.1 y luego de su estimación que se calcula de manera iterativa, los coeficientes son actualizados:

\[\beta = \beta_b+\delta \beta\]

La corrección del bias se calcula utilizanod una regresión lineal con \(N\) predictores independientes \(P_i(x)\) y sus correspondientes coeficientes \(\beta_i\).

\[BC = \sum_{i=1}^{N} \beta_i P_i(x) \]

La corrección del bias luego es aplicada directamente sobre la innovación \(y - H(x)\).

Los predictores utilizados para la corrección asociada al estado de la atmósferea en este sistema corresponden a:

Mientras que la corrección asociada al ángulo de observación utiliza un polinomio de grado N.

El cálculo de CLW y TLR depende de cada sensor y se realiza como parte de las rutinas de GSI.

De acuerdo a Zhu et al. (2019) los coeficientes utilizados en la corrección del bias deben ser entrenados durante un periodo largo de tiempo y muestra que este tiempo es variable para cada canal de cada sensor. Con esto en mente se propone el siguiente experimento para el entrenamiento de los coeficientes.

0.0.1 Diseño experimental

Distintos trabajos abordan el entrenamiento de los coeficientes de corrección de bias en dominios regionales donde la cobertura de observaciones satélitales no es homogenea de manera similar. Por ejemplo Schwartz et al. (2012) utiliza un “spin up” de una semana para el entrenamiento de los coeficientes utilizando un método variacional con ciclos de asimilación cada 6 horas en un dominio regional mientras que Liu et al. (2012) utiliza un método “offline” durante 3 meses para obtener coeficientes de corrección estables. En particular Zhu et al. (2019) propone un periodo de entrenamiento de 9 días utilizando un método variacional determina que es necesario repetir este periodo una vez más para obtener coeficientes estables.

Para obtener coeficientes que permitan la correcta asimilación de observaciones de radianzas durante el periodo 18 UTC 20/11/2018 - 12 UTC 23/11/2018 se realizará asimilación continua entre las 18 UTC 11/11/2018 y las 12 UTC del 20/11/2018.

En este periodo de entrenamiento se observaron distintos fenómenos de tiempo severo asociados a sistemas convectivos de mesoescala y periodos de menor actividad por lo que se espera que la variedad de condiciones atmósfericas favorescan el entrenamiento de coeficientes pese al corto periodo de tiempo (9 días). Para garantizar esto y siguiendo a Zhu et al. (2019) este entrenamiento se realizará dos veces. Para el primer ciclo de asimilación se utilizaron los coeficientes de corrección de bias utilizados por GFS para la misma hora con la intención de acelerar el entramiento en comparación con una inicialización desde “0”. Al finalizar el primer período de asimilación (experimento pre-RAD1) los coeficientes resultantes del último ciclo son utilizados al comenzar el segundo periodo de asimilación en las mismas fechas (pre-RAD2).

El experimento de asimilación se realizará en el dominio que se muestra en la Figura 1 realizando ciclos contínuos a cada hora. El dominio tiene una reticula horizontal de 10 km y el tope del modelo en 50 hPa. En cada ciclo se incluyen las observaciones disponibles en la ventana de asimilación de 1 hora centrada en la hora del análisis incluyendo observaciones convencionales en el prepBUFR, observaciones de estaciones automáticas en el dominio, vientos derivados de satélite y radianzas de sensores en satélites polares listados en la Tabla 1.

Dominio horizontal del modelo

Figure 1: Dominio horizontal del modelo

Table 1: Lista de sensores con la cantidad de canales aceptados por GSI
Sensor Plataforma N° de canales asimilables
airs aqua 118
amsua n15 13
n18 13
aqua 7
n19 13
metop-b 6
atms npp 21
n20 21
cris-fsr npp 100
n20 100
hirs4 metop-a 14
n19 14
iasi metop-a 174
metop-b 174
mhs n18 5
metop-a 5
n19 5
metop-b 5

0.0.1.1 Radianzas

La información de radianzas provienen de los sensores AMSU-A, AIRS, MHS, AIRS, IASI, HIRS todos en satélites de orbita polar. Estas observaciones son parte del sistema de transferencia de datos que utiliza GFS disponibles cada 6 horas pero asimiladas en ciclos horarios utilizando las observaciones correspondientes a una ventana centrada en la hora del análisis.

Se utilizaron observaciones en aire claro (ver métodos para detectar nubles usado por GSI) y se aplicó un control de calidad específico para cada sensor a fin de rechazar observaciones erroneas. A si mismo para tener en cuenta la resolución de cada observación se aplicó un thining de 60 km. Este valor fue definido siguiendo la metodología de otros trabajos en el área (Zhu et al. (2016), Jones et al. (2020), Lin et al. (2017), and Singh et al. (2016)).

En la Figura 2 se muestra la cantidad de observaciones potencialmente asimilables, es decir, que pasan el control de calidad del operador de las observaciones, para cada canal en cada ciclo de asimilación. Dejando de lado IASI que no se ve nada porque tiene más de 600 canales, se puede ver que en todo los casos hay outliers. Esto tiene que ver con que en cada ciclo solo hay disponibles observaciones de algunos sensores y no siempre la pasada del sensor coincide completamente con el dominio. Los momenos donde la cantidad de observaciones aumenta considerablemente coinciden con los momentos donde la pasada del satélite estaba centrada en el dominio. En general hay más observaciones en el rango infrarrojo tanto por las observaciones de los sensores multiespectrales (AIRS e IASI) como por el sensor HIRS que en general para todos los canales aporta mayor cantidad de observaciones.

Cantidad de observaciones potencialmene asimilables en cada ciclo de asimilación para el experimento pre-RAD1

Figure 2: Cantidad de observaciones potencialmene asimilables en cada ciclo de asimilación para el experimento pre-RAD1

En cuanto a la altura de las observaciones (Figura 3), determinada a partir del máximo de la función de peso, muchas observaciones de canales de todos los sensores menos MHS correspoenden a alturas por encima del tope del modelo (50 hPa). Respecto de observaciones en el rango de microondas el sensor ATMS aporta la mayor cantidad de observaciones en niveles medios y altos y lo mismo ocurre con AMSU-A. En cuanto a las observaciones en el infrarrojo AIRS aporta muchas observaciones entre superficie y 200 hPa mientras que los canales de HIRS4 con observaciones en niveles medios y bajos muestran un comportamiento más variane y aporta en general menos observaciones.

TODO: ¿Qué está pasando con el nivel de presión de IASI? Hay una franja negra en 300 hPa que tiene pinta sospechosa.

Nivel de presión asociado a las observaciones potencialmene asimilables en cada ciclo de asimilación para el experimento pre-RAD1

Figure 3: Nivel de presión asociado a las observaciones potencialmene asimilables en cada ciclo de asimilación para el experimento pre-RAD1

0.0.2 Evolución de los coeficientes de corrección.

El entrenamiento de los coeficientes varia sensor a sensor y canal a canal, de acuerdo a lo observado en la bibliografía el tiempo de entrenamiento necesario es variable en cada caso. Al mismo tiempo al ser asimilación regional y horaria podría afectar el tiempo que demora la estabilización de los coeficientes de corrección.

Por ejemplo en la Figura 4 se muestra la evolución de algunos de los coeficientes de canales del AMSU-A en NOAA15 en linea llena para la primera semana de entrenamiento pre-RAD1 y en línea punteada para la segunda semana pre-RAD2. Al igual que Zhu et al. (2019) se observa que cada coeficiente se estabiliza a una velocidad distinta, por ejemplo el coeficiente de emisividad en el canal 8 si bien hacia el final de la semana 2 muestra un comportamiento similar a la semana 1, los valores no son tan cercanos mientras que el TLR tiene un comportamiento muy similar en casi todo el periodo. Lamentablemente Zhu et al. (2019) no muestra otros canales (solo canal 5) pero es esperable que cada uno muestre variaciones.

En general algunos coeficientes para algunos canales no requieren casi entrenamiento mientras que a otros capaz les vendría bien una semanita más.

link a gráficos de todos los sensores y canales

Coeficientes de corrección de bias en cada ciclo de asimilación para los canales 5 a 9 del sensor AMSU-A de NOAA-15.

Figure 4: Coeficientes de corrección de bias en cada ciclo de asimilación para los canales 5 a 9 del sensor AMSU-A de NOAA-15.

Otro ejemplo se observa en la Figura 5 correspondiente a los coeficientes el sensor MHS a bordo del satelite METOP-A. En este caso por ejemplo para el canal 3 se observa como el Offset se mantuvo aprximadamente constante luego de 6 días de entrenamiento mientras que el TLR en el canal 5 continua variando aún luego de la segunda semana.

Coeficientes de corrección de bias en cada ciclo de asimilación para los canales del sensor MHS en METOP-A.

Figure 5: Coeficientes de corrección de bias en cada ciclo de asimilación para los canales del sensor MHS en METOP-A.

0.0.3 Corrección del bias

0.0.3.1 AMSU-A

Si observamos la corrección de bias para cada canal de la familia de sensores AMSU-A se ven diferencias interesantes entre los dos periodos pre-RAD1 y pre-RAD2. Para analizar estas diferencias en la Figura 6 hay que tener en cuenta que solo incluye las observaciones cumplen las condiciones para ser asimiladas y que este BIAS se obtiene de la diferencia entre la innovación con bias corregido y sin corregir. Por lo tanto este es el BIAS corregido por el sistema. Llama la atención que en muchos casos el BIAS durante el periodo pre-RAD2 tiene menor dispersión (menor rango intercuartil) y su mediana está más cerca de cero. Esto podría deberse a que a las observaciones que finalmente pasan el control de calidad requieren una menor corrección o que durante el periodo pre-RAD1 la corrección es sobreestimada.

La evolución de bias a lo largo del tiempo para el sensor sobre NOAA15 (Figura 7) muestra cambios entre los dos periodos en algunos casos del orden de 0.5 grados. En general para los canales 7 y 8 los valores de bias en el periodo pre-RAD2 son algo menores respecto del otro periodo y en los casos donde esto no ocurre as diferencias son muy pequeñas. Esto coincide con lo que se observa en el panel superior derecho de la Figura 7.

Bias de las radianzas de AMSU-A en las distintas plataformas y para los dos periodos.

Figure 6: Bias de las radianzas de AMSU-A en las distintas plataformas y para los dos periodos.

Bias medio de las radianzas de AMSU-A NOAA 15 en las distintas plataformas y para los dos periodos a lo largo del tiempo. El tamaño de los puntos indica la cantidad de observaciones disponibles en cada momento.

Figure 7: Bias medio de las radianzas de AMSU-A NOAA 15 en las distintas plataformas y para los dos periodos a lo largo del tiempo. El tamaño de los puntos indica la cantidad de observaciones disponibles en cada momento.

En la Figura 8 se muestra la contribución de cada predictor a la corrección total de bias, o sea lo que resulta de \(\beta_i P_i(x)\), sin tener en cuenta la corrección asociada al ángulo (TODO: no lo estoy sacando de los diag files originales, capaz debería hacerlo). En primer lugar comparando con el trabajo de Zhu et al. (2019) (Fig. 7) no se ve la variación casi cíclica pero esto muy posiblemente se deba a que en esos experimentos hay observaciones disponibles en cada ciclo de asimilación, lo que no ocurre en los que se presentan acá debido a la asimilación frecuente y al dominio que es más chico. En cuanto a valores, hay mayor similitud aún teniendo en cuenta las diferencias entre los distintos canales.

Analizando específicamente los resultados de los experimentos, se ve que la contribución de cada predictor es distinta para los canales mostrados y además es distinta entre sensores. Por ejemplo Para el canal 5 del AMSU-A NOAA19 el predictor asociado a emisividad es el que más contribuye mientras que para NOAA18 el TLR^2 tambien contribuye. Esto se ve también en canal 7 donde a veces contribuye más el TLR, otras veces el TLR^2.

Evolución temporal del bias asociado a cada predictor para los sensores AMSU-A en las distintas plataformas

Figure 8: Evolución temporal del bias asociado a cada predictor para los sensores AMSU-A en las distintas plataformas

La distribución espacial de la corrección de bias aportada por distintos predictores para un tiempo particular se puede observar en la Figura 9. En este caso el TLR^2 muestra valores comparables con la literatura y en particular un máximo que ocupa gran parte del dominio. El predictor asociado a la emisividad presenta valores distinto de cero solo en la zona de topografía donde el nivel de presión de las observaciones está cerca de superficie.

Contribución a la corrección del bias de dos predictores para el canal 8 de AMSU-A NOAA19 a las 8 UTC del 17/11/2018. Se incluyen observaciones que no pasaron el control de calidad.

Figure 9: Contribución a la corrección del bias de dos predictores para el canal 8 de AMSU-A NOAA19 a las 8 UTC del 17/11/2018. Se incluyen observaciones que no pasaron el control de calidad.

0.0.3.2 MHS

La corrección del bias para el sensor MHS en cada periodo es menos clara que para AMSU-A (Figura 10). Los cambios en pre-RAD2 son muy pequeños y en algunos casos se observan valores muy grandes y mayores a pre-RAD1. Por ejemplo el canal 1 que tiene el máximo cerca de superficie es el que muestra mayor dispersión, esto disminuye mucho para el sensor sobre METOP-A y B (aunque en este ultimo caso la mediana se aleja de cero) pero para NOAA19 la dispersión no cambia particularmente (pero la mediana está más cerca de cero). Esto podría estar indicando que la corrección de bias con los coeficientes entrenados es mayor.

Bias de las radianzas de MHS en las distintas plataformas y para los dos periodos.

Figure 10: Bias de las radianzas de MHS en las distintas plataformas y para los dos periodos.

De la misma manera la evolución temporal del bias para este sensor es muy variable y el periodo pre-RAD2 no muestran cambios marcados (Figura 11. Los valores de bias corregido aumentan y disminuyen alternativamente [TODO: podría deberse a la cantidad de observaciones? eso parece cumplirse en el canal 4 pero se ve lo contrario en el canal 2].

Bias medio de las radianzas de MHS para los dos periodos a lo largo del tiempo.El tamaño de los puntos indican la cantidad de observaciones en cada ciclo.

Figure 11: Bias medio de las radianzas de MHS para los dos periodos a lo largo del tiempo.El tamaño de los puntos indican la cantidad de observaciones en cada ciclo.

En cuanto a la contribución para la corrección del bias por parte de los distintos predictores en MHS (Figura 12) se ven mayores diferencias entre los dos periodos. El canal 1 tiene un pico de 5 K al comienzo del periodo y los canales 4 y 5 tienen algunos valores cercanos a 3 K sobre todo en el periodo pre-RAD2. Para el resto de los canales los valores son similares a lo observado para AMSU-A (Figura 8). El predictor que más contribuye al bias total es distintos para cada canal, por ejemplo en el canal 3 el más importante es el Offset mientras que para el canal 4 es el TLR^2.

Evolución temporal del bias asociado a cada predictor para los sensores MSH en las distintas plataformas

Figure 12: Evolución temporal del bias asociado a cada predictor para los sensores MSH en las distintas plataformas

Para estos sensores en particular llama la atención las diferencias entre los dos periodos, por ejemplo en el canal 4 de METOP-B hay un cambio de magnitud del bias que ocurre alrededor de 16 de noviembre y que se mantiene luego en pre-RAD2 en forma de picos. Para explorar la posibilidad de que este cambio se deba realmente al entrenamiento de los coeficiente y que a partir de esos momentos la corrección del bias esté mejorando, en la Figura 13 se muestra la media de la innovación con bias corregido (línea llena) y bias sin corregir (línea punteada). Se espera que la distribución de la innovación sea aproximadamente gausiana y por lo tanto que la media esté cerca de cero. Una desvíación de cero podría indicar que queda un bias residual luego de la corrección.

En la figura se ven las diferencias bastante marcadas entre la media de la innovación con y sin bias corregido y que en la mayoría de los casos el primer es más cercano a cero, indicando que la corrección del bias va en la dirección correcta. En particular para el canal 4 de METOP-B se obsera que la media de la innovación corregida luego de 16 de noviembre se aleja de cero y esto se mantiene también durante el periodo pre-RAD2. Si bien no es posible identificar ninguna causa esto podía indicar que el coeficiente asociado al TLR^2 cambió afectando negativamente a la corrección del bias y esto se mantuvo a lo largo del tiempo.

Valor medio de la diferencia entre las observaciones y el pronóstico (O - B) a lo largo del tiempo para cada periodo, sensor y canal

Figure 13: Valor medio de la diferencia entre las observaciones y el pronóstico (O - B) a lo largo del tiempo para cada periodo, sensor y canal

La distribución espacial de la corrección del bias asociado a los predictores de TLR^2 y emisividad para un tiempo en particular se muestran en la Figura 14. En este tiempo la pasada del satélite cubrió gran parte del dominio. Respecto de la corrección asociada a la emisividad se observan valores distintos de cero en algunas regiones muy localizadas (sucede parecido en otros tiempos) en las zonas de mayor topografía. Esto podría deberse a que el nivel de presión asociado a esas observaciones se encuentra más cerca de la superficie? El TLR^2 tiene valores pequeños y varian de signo según la región.

Contribución a la corrección del bias de dos predictores para el canal 3 de MHS NOAA19 a las 8 UTC del 17/11/2018. Incluye todas las observaciones de la pasada, aún si no pasan el control de calidad.

Figure 14: Contribución a la corrección del bias de dos predictores para el canal 3 de MHS NOAA19 a las 8 UTC del 17/11/2018. Incluye todas las observaciones de la pasada, aún si no pasan el control de calidad.

0.0.4 Algunos casos de estudio

A continuación se analizan algunos casos de estudio desde el punto de vista de la precipitación acumulada en una hora comparando el análisis con la precipitación estimada por IMERG.

0.0.4.1 2018-11-13 01:00Z

El primer caso corresponder a las 01 UTC del 13 de noviembre y corresponde al 8vo ciclo de asimilación.

En la Figura 15 se muestra arriba la precipitación acumulada para la media del ensamble según cada experimento, abajo a la izquierda la diferencia entre los dos campos y a la derecha la precipitación acumulada según IMERG. En primer lugar es importante notar que los acumulados son mucho menores en los experimentos pero al mismo tiempo hay que tener en cuenta que en la media de los ensables se pierden las características de los miembros con distintas parametrizaciones. En particular algunas parametrizaciones muestran sistemáticamente menor convección afectando la precipitación acumulada.

Las diferencias entre los dos periodos se observan principalmente en toda la franja de precipitación con una disminución en el experimento pre-RAD2 sobre todo en Entre Rios, sur de Uruguay y Rio de la Plata, donde ambos experimentos subestiman la precipitación pero la sobreestiman sobre Santiago del Estero donde no se registró lluvia. Esta sobreestimación disminuye en pre-RAD2 y al mismo tiempo corrige la precipitación sobre Córdoba. Sobre Jujuy y Salta pre-RAD2 muestra mayor precipitación pero ubicada al oeste de las observaciones.

Precipitación acumulada en 60 minutos para cada experimento (arriba), diferencia entre los experimentos (abajo, izquierda) y precipitación acumulada según IMERG (abajo, derecha)

Figure 15: Precipitación acumulada en 60 minutos para cada experimento (arriba), diferencia entre los experimentos (abajo, izquierda) y precipitación acumulada según IMERG (abajo, derecha)

En cuanto a la probabilidad de precipitación acumulada mayor a 5 mm (Figura 16) las diferencias son importantes en la zona del sur de Santa Fe, Entre Rios y Buenos Aires principalmente. En pre-RAD1 la precipitación acumulda parece mucho más extendida en el dominio pero con mayor dispersión entre los miembros del ensable. En pre-RAD2 aparecen valores de probabilidad de entre 30 y 50% sobre la provincia de Córodba y sobre Salta en la región donde se observa precipitación acumulada. Al mismo tiempo sobre Buenos Aires se observa una probabilidad de entre 10 y 20% al sur y mayor al 30% al este que no se observan en la realidad.

Probabilidad de precipitación acumulada > 5 mm (arriba) y > 25 mm (abajo) para cada experimento en el caso 1.

Figure 16: Probabilidad de precipitación acumulada > 5 mm (arriba) y > 25 mm (abajo) para cada experimento en el caso 1.

0.0.4.2 2018-11-18 10:00Z

El segundo caso de estudio corresponde a las 10 UTC del 18 de noviembre casi 7 días luego del comienzo de cada periodo de asimilación continua. En este caso la precipitación se da sobre el norte de argentina y sur de Paraguay y Brasil, con algo de precipitación sobre Salta.

En el campo de diferencia en la Figura 17 se observa que en pre-RAD2 disminuye la precipitación acumulada sobre Salta, Jujuy y Tucumán (esta última no se observa en IMERG) mientras que aumenta sobre Formosa y Chaco donde se dan los mayores acumulados de precipitación. De nuevo, los paneles superiores muestran la media del ensable por lo que muchas características de la ubicación de la precipitación acumulada se pierden y hay mucha variabilidad entre los miembros.

Precipitación acumulada en 10 minutos para cada experimento (arriba), diferencia entre los experimentos (abajo, izquierda) y precipitación acumulada en 30 minutos según IMERG (abajo, derecha)

Figure 17: Precipitación acumulada en 10 minutos para cada experimento (arriba), diferencia entre los experimentos (abajo, izquierda) y precipitación acumulada en 30 minutos según IMERG (abajo, derecha)

Esto puede verse en la Figura 18 que muestra la probabilidad de precipitación acumulada mayor a 5 mm. En pre-RAD2 se ve mayor probabilidad de precipitación por encima de 50% extendida sobre gran parte de la provincia de Formosa mientras que la precipitación observada según IMERG se localiza más hacia el este. La precipitación que se observa sobre Misiones en los experimentos no debe ser considerada con mucha confianza ya que se encuentran en el borde del diminio de simulación.

Probabilidad de precipitación acumulada > 5 mm (arriba) y > 25 mm (abajo) para cada experimento correspondiente al caso 2

Figure 18: Probabilidad de precipitación acumulada > 5 mm (arriba) y > 25 mm (abajo) para cada experimento correspondiente al caso 2

0.0.5 Colophon

This report was generated on 2020-10-21 21:10:33 using the following computational environment and dependencies:

#> ─ Session info ───────────────────────────────────────────────────────────────
#>  setting  value                         
#>  version  R version 3.6.3 (2020-02-29)  
#>  os       Debian GNU/Linux 9 (stretch)  
#>  system   x86_64, linux-gnu             
#>  ui       X11                           
#>  language (EN)                          
#>  collate  en_US.UTF-8                   
#>  ctype    en_US.UTF-8                   
#>  tz       America/Argentina/Buenos_Aires
#>  date     2020-10-21                    
#> 
#> ─ Packages ───────────────────────────────────────────────────────────────────
#>  package            * version    date       lib
#>  assertthat           0.2.1      2019-03-21 [1]
#>  backports            1.1.10     2020-09-15 [1]
#>  bitops               1.0-6      2013-08-17 [1]
#>  bookdown             0.20       2020-06-23 [1]
#>  broom                0.5.6      2020-04-20 [1]
#>  callr                3.4.4      2020-09-07 [1]
#>  cellranger           1.1.0      2016-07-27 [2]
#>  checkmate            2.0.0      2020-02-06 [1]
#>  class                7.3-15     2019-01-01 [4]
#>  classInt             0.4-3      2020-04-07 [1]
#>  cli                  2.0.2      2020-02-28 [1]
#>  colorspace           1.4-1      2019-03-18 [1]
#>  crayon               1.3.4      2017-09-16 [2]
#>  data.table         * 1.12.8     2019-12-09 [1]
#>  DBI                  1.1.0      2019-12-15 [1]
#>  dbplyr               1.4.3      2020-04-19 [1]
#>  desc                 1.2.0      2018-05-01 [2]
#>  devtools             2.3.2      2020-09-18 [1]
#>  digest               0.6.25     2020-02-23 [1]
#>  dplyr              * 1.0.0      2020-05-29 [1]
#>  e1071                1.7-3      2019-11-26 [1]
#>  ellipsis             0.3.1      2020-05-15 [1]
#>  evaluate             0.14       2019-05-28 [1]
#>  fansi                0.4.1      2020-01-08 [1]
#>  farver               2.0.3      2020-01-16 [1]
#>  forcats            * 0.5.0      2020-03-01 [1]
#>  fs                   1.5.0      2020-07-31 [1]
#>  generics             0.0.2      2018-11-29 [1]
#>  ggplot2            * 3.3.1      2020-05-28 [1]
#>  glue                 1.4.2      2020-08-27 [1]
#>  gtable               0.3.0      2019-03-25 [1]
#>  haven                2.2.0      2019-11-08 [1]
#>  here               * 0.1        2017-05-28 [1]
#>  highr                0.8        2019-03-20 [1]
#>  hms                  0.5.3      2020-01-08 [1]
#>  htmltools            0.5.0      2020-06-16 [1]
#>  httr                 1.4.2      2020-07-20 [1]
#>  isoband              0.2.1      2020-04-12 [1]
#>  jsonlite             1.7.1      2020-09-07 [1]
#>  kableExtra         * 1.2.1      2020-08-27 [1]
#>  KernSmooth           2.23-16    2019-10-15 [4]
#>  knitr              * 1.30       2020-09-22 [1]
#>  labeling             0.3        2014-08-23 [2]
#>  lattice              0.20-38    2018-11-04 [4]
#>  lifecycle            0.2.0      2020-03-06 [1]
#>  lubridate          * 1.7.9      2020-06-08 [1]
#>  magrittr             1.5        2014-11-22 [2]
#>  memoise              1.1.0      2017-04-21 [2]
#>  mesoda             * 0.0.0.9000 2020-10-08 [1]
#>  metR               * 0.8.9000   2020-10-16 [1]
#>  modelr               0.1.7      2020-04-30 [1]
#>  munsell              0.5.0      2018-06-12 [2]
#>  ncdf4                1.17       2019-10-23 [1]
#>  nlme                 3.1-143    2019-12-10 [4]
#>  patchwork          * 1.0.0      2019-12-01 [1]
#>  pillar               1.4.6      2020-07-10 [1]
#>  pkgbuild             1.1.0      2020-07-13 [1]
#>  pkgconfig            2.0.3      2019-09-22 [1]
#>  pkgload              1.1.0      2020-05-29 [1]
#>  plyr                 1.8.6      2020-03-03 [1]
#>  prettyunits          1.1.1      2020-01-24 [1]
#>  processx             3.4.4      2020-09-03 [1]
#>  proj4                1.0-10     2020-03-02 [1]
#>  ps                   1.3.4      2020-08-11 [1]
#>  purrr              * 0.3.4      2020-04-17 [1]
#>  R6                   2.4.1      2019-11-12 [1]
#>  RColorBrewer         1.1-2      2014-12-07 [2]
#>  Rcpp                 1.0.5      2020-07-06 [1]
#>  RCurl                1.98-1.2   2020-04-18 [1]
#>  readr              * 1.3.1      2018-12-21 [1]
#>  readxl               1.3.1      2019-03-13 [1]
#>  Redmonder            0.2.0      2017-01-04 [1]
#>  remotes              2.2.0      2020-07-21 [1]
#>  reprex               0.3.0      2019-05-16 [1]
#>  rgeos                0.5-2      2019-10-03 [1]
#>  rlang                0.4.7      2020-07-09 [1]
#>  rmarkdown            2.4        2020-09-30 [1]
#>  rnaturalearth        0.1.0      2017-03-21 [1]
#>  rnaturalearthhires   0.2.0      2020-01-23 [1]
#>  rprojroot            1.3-2      2018-01-03 [2]
#>  rstudioapi           0.11       2020-02-07 [1]
#>  rvest                0.3.5      2019-11-08 [1]
#>  scales               1.1.1      2020-05-11 [1]
#>  selectr              0.4-2      2019-11-20 [1]
#>  sessioninfo          1.1.1      2018-11-05 [1]
#>  sf                   0.9-2      2020-04-14 [1]
#>  sp                   1.4-1      2020-02-28 [1]
#>  stringi              1.5.3      2020-09-09 [1]
#>  stringr            * 1.4.0      2019-02-10 [1]
#>  testthat             2.3.2      2020-03-02 [1]
#>  tibble             * 3.0.3      2020-07-10 [1]
#>  tidyr              * 1.1.0      2020-05-20 [1]
#>  tidyselect           1.1.0      2020-05-11 [1]
#>  tidyverse          * 1.3.0      2019-11-21 [1]
#>  unglue             * 0.0.1      2019-10-07 [1]
#>  units                0.6-6      2020-03-16 [1]
#>  usethis              1.6.3      2020-09-17 [1]
#>  vctrs                0.3.4      2020-08-29 [1]
#>  viridisLite          0.3.0      2018-02-01 [2]
#>  webshot              0.5.2      2019-11-22 [1]
#>  withr                2.3.0      2020-09-22 [1]
#>  xfun                 0.18       2020-09-29 [1]
#>  xml2                 1.3.2      2020-04-23 [1]
#>  yaml                 2.2.1      2020-02-01 [1]
#>  source                                      
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.3.3)                              
#>  local                                       
#>  Github (eliocamp/metR@385c262)              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  Github (ropensci/rnaturalearthhires@2ed7a93)
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.3.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.3)                              
#>  CRAN (R 3.6.2)                              
#>  CRAN (R 3.6.2)                              
#> 
#> [1] /home/paola.corrales/R/x86_64-pc-linux-gnu-library/3.6
#> [2] /usr/local/lib/R/site-library
#> [3] /usr/lib/R/site-library
#> [4] /usr/lib/R/library

The current Git commit details are:

#> Local:    main /home/paola.corrales/mesoda
#> Remote:   main @ origin (git@github.com:paocorrales/mesoda.git)
#> Head:     [b80f298] 2020-10-20: Bias residual

Jones, Thomas A., Patrick Skinner, Nusrat Yussouf, Kent Knopfmeier, Anthony Reinhart, Xuguang Wang, Kristopher Bedka, William Smith, and Rabindra Palikonda. 2020. “Assimilation of GOES-16 Radiances and Retrievals into the Warn-on-Forecast System.” Monthly Weather Review 148 (5): 1829–59. https://doi.org/10.1175/MWR-D-19-0379.1.

Lin, Haidao, Stephen S. Weygandt, Agnes H. N. Lim, Ming Hu, John M. Brown, and Stanley G. Benjamin. 2017. “Radiance Preprocessing for Assimilation in the Hourly Updating Rapid Refresh Mesoscale Model: A Study Using AIRS Data.” Weather and Forecasting 32 (5): 1781–1800. https://doi.org/10.1175/WAF-D-17-0028.1.

Liu, Zhiquan, Craig S. Schwartz, Chris Snyder, and So-Young Ha. 2012. “Impact of Assimilating AMSU-A Radiances on Forecasts of 2008 Atlantic Tropical Cyclones Initialized with a Limited-Area Ensemble Kalman Filter.” Monthly Weather Review 140 (12): 4017–34. https://doi.org/10.1175/MWR-D-12-00083.1.

Miyoshi, Takemasa, Yoshiaki Sato, and Takashi Kadowaki. 2010. “Ensemble Kalman Filter and 4D-Var Intercomparison with the Japanese Operational Global Analysis and Prediction System.” Monthly Weather Review 138 (7): 2846–66. https://doi.org/10.1175/2010MWR3209.1.

Schwartz, Craig S., Zhiquan Liu, Yongsheng Chen, and Xiang-Yu Huang. 2012. “Impact of Assimilating Microwave Radiances with a Limited-Area Ensemble Data Assimilation System on Forecasts of Typhoon Morakot.” Weather and Forecasting 27 (2): 424–37. https://doi.org/10.1175/WAF-D-11-00033.1.

Singh, Randhir, Satya P. Ojha, C. M. Kishtawal, P. K. Pal, and A. S. Kiran Kumar. 2016. “Impact of the Assimilation of INSAT-3D Radiances on Short-Range Weather Forecasts: Assimilation of INSAT-3D Radiances.” Quarterly Journal of the Royal Meteorological Society 142 (694): 120–31. https://doi.org/10.1002/qj.2636.

Zhu, Kefeng, Ming Xue, Yujie Pan, Ming Hu, Stanley G. Benjamin, Stephen S. Weygandt, and Haidao Lin. 2019. “The Impact of Satellite Radiance Data Assimilation Within a Frequently Updated Regional Forecast System Using a GSI-Based Ensemble Kalman Filter.” Advances in Atmospheric Sciences 36 (12): 1308–26. https://doi.org/10.1007/s00376-019-9011-3.

Zhu, Yanqiu, Emily Liu, Rahul Mahajan, Catherine Thomas, David Groff, Paul Van Delst, Andrew Collard, Daryl Kleist, Russ Treadon, and John C. Derber. 2016. “All-Sky Microwave Radiance Assimilation in NCEP’s GSI Analysis System.” Monthly Weather Review 144 (12): 4709–35. https://doi.org/10.1175/MWR-D-15-0445.1.